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ABSTRACT 
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CN . We studied the relation between the ratio of rotational velocity to velocity 

dispersion and the metallicity (< v rot > /ovmetallicity relation) of globular 
\ cluster systems (GCS) of disk galaxies by comparing the relation predicted from 

CO \ simple chemo-dynamical models for the formation and evolution of disk galaxies 

i | with the observed kinematical and chemical properties of their GCSs. 

We conclude that proto disk galaxies underwent a slow initial collapse that 
\ was followed by a rapid contraction and derive that the ratio of the initial collapse 

time scale to the active star formation time scale is ~ 6 for our Galaxy and ~ 15 
for M31. 

The fundamental formation process of disk galaxies was simulated based on 
simple chemo-dynamical models assuming the conservation of their angular mo- 

Q-V mentum. We suggest that there is a typical universal pattern in the < v rot > /o~ v - 

metallicity relation of the GCS of disk galaxies. This picture is supported by the 
observed properties of GCSs in the Galaxy and in M31. This relation would 
deviate from the universal pattern, however, if large-scale merging events took 
major role in chemo-dynamical evolution of galaxies and will reflect the epoch of 
such merging events. We discuss the properties of the GCS of M81 and suggest 

c3 ' the presence of past major merging event. 

Subject headings: galaxies: formation, galaxies: evolution, galaxies: star clusters, 
globular clusters 



1. Introduction 

The importance of investigations of globular clusters (GCs) has been reassessd in recent 
years (Ashman & Zepf 1998, Kissler-Patig 2000), particularly since the globular cluster 
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systems (GCS) of disk galaxies can be used to study the dynamical formation history of 
their host galaxies. 

The rapid dynamic formation scenario of our Galaxy was first proposed by Eggen, 
Lynden-Bell & Sandage (1962; hereafter ELS) based on apparent correlations among the 
kinematic properties and metallicity indicators of the halo and disk stars in the solar neigh- 
borhood. 

Previous studies interpreted the properties of the GCS in relation to the formation of 
our Galaxy. Maline, Hartmann & Mathews (1991) treated the GCs of the Galaxy as a tracer 
of its chemical evolution and concluded that the Galaxy was formed by an inhomogeneous 
collapse. 

Kinman (1959) demonstrated that the GCS of our Galaxy as a whole has a small 
net rotation and Frenk & White (1980) suggested that the difference in rotational velocity 
between the F-type (metal-poor) and G-type (metal-rich) GC groups is only marginally 
significant in our Galaxy. Huchra et al. (1982) also discovered a rotating system of M31 
GCs. Zinn (1993) pointed out that inner "old halo" clusters have rapid rotation, small 
velocity dispersion, and flattened orbits and suggested that these properties are evidence of 
the dissipational collapse formation of our Galaxy. 

Minniti (1996) studied the kinematics of GCs and stars in the Galactic bulge and con- 
cluded that the metal-rich GCs and stars in the Galactic bulge have similar kinematic prop- 
erties to the inner "old halo" population (cf. Zinn 1993) and that the Galactic bulge was 
formed by dissipational collapse. 

Recently, Chiba & Beers (2000) discussed a new scenario for the formation of our Galaxy 
from studying the proper motions of metal-poor halo stars in the solar neighborhood obtained 
by HIPPARCOS. 

In studying the global dynamic evolution of galaxies, it is advantageous to use the 
kinematic properties of the GCSs of those galaxies, because most GCs should have been 
formed before, during, or immediately after the collapse of their host galaxies (Fall & Rees 
1985). Those GCs, therefore, can give us information on the dynamic evolution that occurred 
during those epochs. 

We discuss the dynamic formation scenario of disk galaxies using the kinematic property 
of the GCS, especially the relation between metallicity and the ratio of the rotational velocity 
to the velocity dispersion. 

We construct simple models to simulate the evolution of the chemo-dynamic properties 
of the GCSs of disk galaxies in Section 2 and compare the observed properties of the GCSs of 
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our Galaxy and M31 with those derived from these simple models in Section 3. In Section 4, 
we discuss the importance of observing GCSs in SO and Sa type galaxies. Finally, in Section 
5, we try to recompile the GCS data of M81 and propose that the chemo-dynamic evolution 
of the M81 system may differ from those of the Galaxy and M31. We then conclude that we 
must investigate the detailed kinematics of the M81 GCS, because the quality and quantity 
of the data on M81 GCs is poor. 



2. Simple Models for the Dynamic Evolution of GCS 

2.1. Model Construction 

In this section, we investigate the time variation of the rotational velocity of the col- 
lapsing protogalactic gas of the disk galaxy from which GCs were formed, and search for a 
relation between the metallicity Z and the normalized specific rotation < v° b l(Z) > ja v iZ) 
of GCS, where < v°^ s (Z) > and a v (Z) are the mean rotational velocity and the velocity 
dispersion of the GCS, respectively. 

In order to make the dynamic models of disk galaxy formation treatable, we assume 
that the protogalaxies underwent homologous dynamic contraction along a family of self- 
consistent density-potential pair solutions of the Poisson equation by conserving its angular 
momentum. 

The axisymmetric models adopted in this paper are those derived by Miyamoto & Nagai 
(1975) and their density p and the gravitational potential \1/ are given by 

p(r, z; ait), bit)) — 1 — —, (1) 

*(r, z- a(t), b(t)) = - - GMga l = , (2) 

yV 2 + (a + Vz 2 + b 2 ) 2 

where M gai is the total mass of the galactic disk and the parameters ait) and bit) are 
the time-dependent scale lengths along the r and z directions, respectively. Solutions (1) and 
(2) for the density-potential pair satisfy the Poisson equation given by 



V 2 ^ = 4ttGp . (3) 
For simplicity, we assume that b/a = 0.2 at all times, which is acceptable for typical 
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disk galaxies. The contraction profile of the protogalaxies is then expressed by the variation 
of the scale length a(t) as a function of time t. 

Furthermore, we assume that the total angular momentum of the protogalaxies is 
conserved during the collapse phase and its direction is perpendicular to the galactic disk. 
The expression for the total angular momentum is given by 



where i^(r, t) and E(r, t) are the rotation velocity parallel to the galactic disk and the 
projected surface density defined by 



where p(r, z, t) is the volume density of the galaxy. For practical calculation purposes, 
we adopt the radius R = 4a and the half thickness H = 46. It has been shown that ~ lCT 6 p c 
at these radius and half thickness defined by R and H, where the p c is the central mass 
density of a disk galaxy. 

Since a flat rotation curve is a typical property of disk galaxies, we assume that 



This requires that when the density profile varies with the time-dependent scale radius 
a(t), the angular momentum distribution is adjusted to make up the flat v^, whereas the 
total angular momentum of the galaxy is conserved. 

The flat rotation velocity is then given by 




(4) 




(5) 



< v^,(r,t) >= V f i at {t) 



(6) 



Vflat{t) 



(7) 



where the moment of inertia Ir(£) is defined by 




(8) 
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For simplification, in this paper, we assume that the spatial distribution, angular mo- 
mentum, total energy, and flat property of the rotation curves of these GCs have been frozen 
and conserved since the time of GCS formation. Furthermore, we assume that once GCs 
form in the collapsing protogalactic gas, they are disconnected from the collapsing gas sys- 
tem. If the kinematical energy of GCS were not conserved, the < v° h Q l(Z) > should be larger 
than that of conserved case for metal-poor clusters. However the a v also should be larger, 
so the effect of the dissipation for < v° b l(Z) > /a v (Z) can be negligible. 

Under these assumptions, the mean rotational velocity < v^(r, Z) > of GCs observed 
at present (t = t\) would be equivalent to the mean rotational velocity Vfi a t(t) during the 
formation of those GCs. The age-metallicity relation that translates t to Z is discussed in 
Section 2.3. 

Other velocity components t>' 1 (r,.Z') and u* 1 ^, Z) depend on < v^(r, Z) > and the 
present gravitational potential \I/(r, z; a(ti), b(ti)). If the disk galaxies are well virialized, 
those velocity components are given by 



.obs 

'rot 



v t /(r,Z)=v t /(r,Z) = ^ 1 -^ f) . (9) 

We then derive the formulation for the observational, normalized specific rotation < 
(Z) > I cr v (Z) of an disk galaxy from those three velocity components and the inclination % 
of the disk projected onto the celestial sphere. By contrast, the observed < v° b s t (Z) > /a v (Z) 
of our Galaxy is derived by assuming that the observer is at the solar radius (see Appendix). 

The parameters adopted in this paper are: the scale length a\ of a typical disk galaxy 
at present is a\ ~ 6.25kpc, hence R\ = 25kpc (e.g. Helmi & White 1999, Patsis & Grosb0l 
1996), the initial (t = t ) disk radius of a galaxy is R = 3Ri (e.g., Rees & Ostriker 1977), 
and the present flat rotation velocity is v rot (ti) ~ 250km s^ 1 (e.g., Rubin & Ford 1970, 
Burton & Gordon 1978). 



2.2. Collapse Profile of Scale Radius a 

It is not straightforward to solve the time variation of the scale radius a(t) of the 
gas self-consistently. Instead, we adopt three models for describing the collapse profiles of 
protogalaxies for consideration. 

(1) Exponential collapse: In this model, we assume that the scale size of the protogalaxy 
decreases exponentially with time, and that a(t) can be written in the form 



6 



a(t) = a exp 




(10) 



where t co \ is the timescale of collapse and ao is the initial size of the protogalaxy. Of 
the three models, this represents a moderate collapse. 

(2) Delayed collapse: In this model, the size of the protogalaxy decreases gradually 
with time initially and then, after some epoch, it decreases more rapidly with time. For our 
Galaxy, based on Zinn(1993), Chiba & Beers (2000), this scenario assumes that the halo 
was built up by numerous fragments falling into the Galaxy (e.g., Searle & Zinn 1978) and 
that at some time the halo collapsed rapidly (e.g., ELS). We assume that the slow collapse 
process was dominant until tt ra ns and that the rapid collapse process became dominant after 
t trans , namely 



where t co i t i and t co / i2 are the long and short time scales (i.e., t co i t \ > t trans > t co i t2 ), 
characterizing the two-stage collapse taking place before and after the transition epoch t trans . 

(3) Early collapse: In this model, the protogalaxy experiences a rapid initial collapse, 
followed by a slow collapse. The equation describing this collapse profile is the same as for 
the delayed collapse model; however, the time scale of the first collapse, t C oZ,n is smaller than 
that of the second collapse, t co i t 2- Although this model is not supported by theoretical or 
observational evidence, we examine this model as an extreme case, in contrast to the delayed 
collapse model. 

The collapse profiles describing the time variation of a(t) for the protogalaxy in these 
three collapse models are shown in Figure 1. 



This subsection discusses the temporal evolution of metallicity, that is, the so-called 
age metallicity relation (AMR) of disk galaxies under the assumption that the metallicity 
of a GC traces the chemical evolution of the collapsing protogalaxy. The age of a galaxy is 
defined as the elapsed time from the start of the initial collapse of the galaxy. The initial 
mass function and rate of star formation are probably dependent on the density and pressure 
of the collapsing gas. Therefore, the AMR might well be modified according to the phase of 



a(t) 
a(t) 




(11) 



2.3. Age Metallicity Relation 
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star formation. For simplification, however, we assume that the AMR is independent of the 
dynamic effect. 

We adopt a simple closed-box model for the chemical evolution of disk galaxies. This 
simple closed-box model is based on Hartwick's (1976) concept, introduced as a simple model 
for the chemical evolution of the Galactic halo. We assume that the GCs were formed from 
the collapsing gas and that the GCs conserve their initial chemical composition. Under this 
assumption, the chemical abundance in GCs should reflect the chemical abundance in the 
collapsing halo gas during the formation of those GCs. 

This simple closed-box model adopts the Schmidt Law ip s f(t) oc f£(t) for the star 
formation rate, where f g (t) is the gas fraction of the galaxy. 

The AMR of a disk galaxy is obtained by coupling the basic equations for chemical 
evolution (Tinsley 1980) with the simple closed-box model. It is given by 

log(Z/Z ) = log(j//Z ) + log(t/f a/ ) (n = 1) , (12) 

or 

log(Z/Z ) = log(y/Z ) + log{log(l + t/t sf )} (n = 2) , (13) 

where Z Q is the metallicity of the solar neighborhood, t is the age of the galaxy since 
the onset of the initial collapse, t s f is the time scale of the active star formation period, and 
y is the yield of newly created elements. These equations are quoted from Arimoto, Yoshii 
& Takahara (1992). For our Galaxy, the observation implies y ~ 0.1Z Q in the halo. In this 
work, we assume that nearby disk galaxies are the same age as the Galaxy, so they have the 
same star formation rate as the Galaxy. 

2.4. Implication of Simple Model 

Figure 1 shows the collapse profiles of the scale radius a(t) for the three models studied 
in this paper. Figure 2 shows the evolution of the normalized specific rotation < v° b f(Z) > 
/a v (Z) of the disk galaxy as a function of metallicity log Z/Z Q for each of the three collapse 
models for n — 1. The exponential and delayed collapse models both show rather an abrupt 
increase in < v° b J(Z) > /a v (Z) at logZ/Z Q ~ -1 for t col /t sf = t trans /t sf = 1 and at 
log Z/Z Q ~ —0.0 for tcoi/tsf = Urans/tsf = 10. The epoch of this abrupt increase in < 
v° b Q l(Z) > /<j v (Z) depends directly on the value of t s f. For the early collapse model, < 
v° b Q l(Z) > /a v {Z) increases gradually in contrast to the exponential and delayed collapse 
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models. Thus, we can discriminate between the early collapse model and the other two 
models in terms of the increase in < v° b f(Z) > /a v (Z). 

Figure 3 shows a relation similar to that of Figure 2, but for n — 2. Again, we see an 
abrupt increase in < v° b l(Z) > /a v (Z) at some log Z/Z Q for the exponential and delayed 
collapse models. In these two models, the abrupt increase in < v° b a l(Z) > /a v (Z) is a 
universal property in contrast to the early collapse model. Therefore, we can easily see 
whether the early collapse model adequately describes the dynamic collapse history of disk 
galaxies, by examining the observed < v° h l(Z) > /a v (Z) of the GCSs in disk galaxies. 

From Figures 2 and 3, we conclude that when n and t co i/t s f or t trans /t s f are fixed in a 
model, the profiles of the increase in < v° b f(Z) > /a v (Z) are very similar. In summary, 

(1) The evolution of the normalized specific rotation as a function of the metallicity 
of GCS, described by the relation < v° b l(Z) > /a v (Z) can be used as a diagnostic tool to 
identify the collapse history of the halos of disk galaxies. 

(2) When the Schmidt Law Index n is fixed, the Z-dependence of the ratio of < v ° h l{Z) > 
/a v (Z) is sensitive to the ratio t co i/t s f, but is insensitive to the individual value of t co i or t s f 
for a common ratio of t co i/t s f. 

3. The Observational Relation between < v° b Q l > jo v and [Fe/H] 

3.1. Our Galaxy 

We now compare the observed kinematic properties of GCSs with those predicted from 
the simple models discussed in Section 2. The data of the GCS of our Galaxy, which are 
used to derive the < v rot > / cvmetallicity relation presented in this paper, are based on the 
most recent catalogue of Harris (1996). In this catalogue, only the heliocentric velocity Vh 
is available. First, we correct Vh to v ra d — Vh + t>© sin /cos b, where v Q is the velocity of the 
LSR about the galactic center, and I and b are Galactic longitude and latitude, respectively. 
Therefore, we regard the mean radial velocity as the rotational velocity. In order to derive it, 
we divide the Galactic disk into two regions: one includes longitudes with 0° < I < 180° and 
the other includes longitudes with 180° < I < 360°. Then, we calculate the mean of radial 
velocity for each of the regions, < v^ ad > and < v~ ad > and regard < v^ ad > — < v~ ad > as 

< v° b Q l >. We use [Fe/H] as an indicator of the metallicity Z. 

We adopt the "running mean" method to evaluate < v° b l(Z) > /a v (Z) from a limited 
number of available samples, and hence are interested only in global properties of the derived 

< v°H{Z) > /a v (Z) . 
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The general properties of the derived < v rot > /ovmetallicity relation of the GCS of 
our Galaxy (Figure 4) are as follows: 

(1) The GCs having [Fe/H] < —1.0 show only a modest net rotation < v rot > /a ~ 0.5. 

(2) The rotational velocity of GCs increases rather abruptly from [Fe/H] ~ —1.0 to 
[Fe/H] ~ -0.7. 

In our Galaxy the metallicity distribution of the GCS shows the presence of two distinct 
components (Zinn 1985). This bimodality suggests that two different formation processes 
were involved. It appears that one of the populations was mixed with the other at [Fe/H]> 
— 1.0. Therefore, the assumption that the GCS was formed from the collapsing proto-Galaxy 
might not hold in that range of metallicity. 

This paper examines the chemical evolution of GCs that formed during the halo col- 
lapse phase with various collapse profiles, but it does not consider the deviation of the 
age-metallicity relation of the disk population of GCs from that of the halo population. Our 
simple model reproduces the general increasing trend of < v° b l{Z) > /cr v (Z) with increasing 
[Fe/H], but the low values of < vfj(Z) > /a v (Z) were reproduced at -1.8 <[Fe/H]< -1.0. 
This might be due to a major merging episode that took place during this epoch, since 
major merger events can increase a v temporally and decrease net < v° b J(Z) > to lower 
< vfj(Z) > /a v (Z) . 



3.2. M31 

M31, the nearest spiral galaxy, has a similar mass (1.8 x 1O 11 M ; Rubin & Frod 1970) 
and morphological type (Sb; de Vaucouleurs et al. 1991) to our Galaxy, and is a natural test 
bed to apply our scheme for studying chemo-dynamic evolution using its GCS. 

We recompiled the most recent GCS catalogue of M31 (Barmby et al. 2000) to derive 
< v rot > /cTy-metallicity relation. We regard the mean radial velocity as the rotational 
velocity. In order to derive it, we divide the M31 disk into two regions: one has a direction 
along with major-axis of M31 disk from center of M31 and the other has the other direction. 
Then, we calculate the mean of radial velocity for each of the regions, < v^ ad > and < v~ ad > 
and regard < v+ ad > - < v~ ad > as < v° b J >. 

By adopting the inclination of M31, 77° (e.g., Rubin & Ford 1970), we derived Figure 
5 for M31, in the same manner as Figure 4. The overall property of < v ro t > /ovmetallicity 
relation for M31 is similar to that for our Galaxy. By comparing the observed < v rot > jo v - 
metallicity relations of the GCSs of our Galaxy and M31 with the results of the simple model 
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discussed in Section 2, we arrived at the following two conclusions: 

(1) The contraction of the protogalaxy of M31 followed a collapse profile similar to all 
collapse models. 

(2) The ratio of the collapse time scale to the star formation time scale should have 
been t co i/t s f ~ 3 — 11 for the Galaxy, while it should be t co i/t s f > 7 for M31, when the power 
index of the Schmidt Law for star formation is n — 1. For n — 2, it should be t co i/t s f > 600 
for the Galaxy and t co i/t s f > 1.0 x 10 4 for M31. However, we notice that t co i/t s f > 100 may 
be unusual value for that ratio. 

The range of t co i/t s f is derived considering uncertainties of our model and the observed 
< v r ot > /er„-metallicity relation, i.e. if a model line with a t co i/t s f is covered error bars of 
metal-rich GCs' data points, this t co i/t s f is accepted. We regard the range of —1.0 <[Fe/H]< 
—0.3 and —0.6 <[Fe/H] as metal-rich GCs' data points for our Galaxy and M31, respectively. 
For our Galaxy, we leave the point of [Fe/H]= —0.2 out of consideration, because it should 
be unexpected point for our model. 

4. Effects of Late Mergers on < v rot > /ovmetallicity relation 

Recently, the hierarchical clustering model has been discussed as a standard galactic 
formation model (e.g., Katz & Gunn 1991). Merger events probably played a considerable 
role in triggering GCs formation (Ashman & Zepf 1992). There is evidence that some 
elliptical galaxies have undergone mergers (e.g., Zepf & Ashman 1993). Zepf & Ashman 
(1993) predicted that the GCSs of some elliptical galaxies have two peaks in their metallicity 
distribution. Our simple model assumes that galaxies were formed through the monolithic 
collapse of primordial gas. Our dynamical model is compatible with a hierarchical clustering 
model, if the clustering gas clumps are numerous and have small scale. However, if the 
clustering of gas clumps involves a major merger event, our model may not adequately 
describe the dynamic and chemical evolution of such a system. 

By contrast, for those dissipationless merger events that did not trigger the formation 
of GCs, only the kinematics of GCs were affected. In fact, Kissler-Patig & Gebhardt (1998) 
found evidence for dissipationless merger events in the globular clusters kinematics of M87. 
Some SO and Sa galaxies might also have undergone dissipationless merger events (e.g., 
Jore, Broeils & Haynes 1996, Mihos, Walker & Hernquist 1995 ). The occurrence of similar 
dissipationless merger events is conceivable in disk galaxies. If dissipationless merger events 
played a significant role in disk galaxies, < v rot > / ovmetallicity relation cannot be predicted 
using the simple models introduced in Section 2. In other words, if the formation of the GCS 
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is affected by dissipationless merging events, we will find that its kinematical properties are 
more complicated and will not show the rather simple relations observed in our Galaxy and 
M31 (cf. Figures 4 and 5). 

We need to investigate the relation between the kinematics and metallicity for the GCSs 
of early-type disk galaxies to see the extent of the role of dissipationless merging events in 
the formation of GCSs. 

5. A Preliminary Assessment of the < v rot > /o\,-metallicity relation for M81 

Elaborate observational studies of the GCSs of disk galaxies have been carried out 
during the last decade. After the Galaxy and M31, the GCS of M81 is one of the most 
extensively observed (Georgiev et al. 1992, 1993, Perelmuter & Racine 1995, Perelmuter, 
Brodie & Huchra 1995). Perelmuter, Brodie & Huchra (1995) published kinematical data 
of the GCs of M81. Even so, at present only 25 M81 GCs are available for the study of 
their chemo-dynamical properties. Figure 6 shows the < v ro t > /ovmetallicity relation 
for M81, indicating a increasing < v° b l(Z) > /a v (Z) at all [Fe/H] similar to our Galactic 
and M31 GCS. M81 likely underwent a significant interaction with M82 and NGC 3077, 
the massive member galaxies of the M81 group. The evidence for such an interaction is 
the neutral hydrogen stream between those galaxies that has been confirmed by 21cm HI 
emission studies (van der Hulst 1979, Yun et al. 1993). Therefore, the kinematics of the 
M81 GCS might have been modified by tidal effects more conspicuously than in the Galaxy 
or M31. 

At present, our study is limited by the paucity of available data on the GCS of M81. 
Since the expected number of GCs in M81 is more than 150, we are planning an extensive 
study of the properties of the GCS of M81 to constrain the scenario of M81 formation. 

6. Conclusion 

We made a series of chemo-dynamic model calculations to predict the relation for GCSs 
under the assumption that GCSs are formed before or during the collapse of protogalactic 
gas clouds. 

Three collapse profiles, e.g., exponential collapse, delayed collapse, and early collapse, 
are assumed in our simple models. The power index n of the Schmidt Law for the star 
formation process, the dynamic collapse time scale t co i, and the star formation time scale t s f 
are the free parameters chosen in this study. 
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We analyzed the < v rot > /<7„-metallicity relations of our Galaxy (Figure 4) and M31 
(Figure 5). The observed relations for these two disk galaxies of similar mass and mor- 
phological type are similar. They exhibit a rather sudden increase in < v° h Q l(Z) > /a v (Z) 
at [Fe/H]~ —1.0 for the Galaxy and at [Fe/H]~ —0.6 for M31. This is consistent with 
chemo-dynamic model predictions of the formation of the GCSs of these galaxies. However 
we have to notice that it is difficult to choose the parameter set (t co i t i,t co i t 2) for fitting the 
model with data points for early collapse model. It also suggests that the GCSs of these 
galaxies were formed in a dynamic collapse, where the collapse time scale is as large as, or 
larger than, the star formation time scale, if the star formation process is proportional to 
the protogalactic gas density, namely n = 1. The collapse time scale should be at least an 
order of magnitude larger than the star formation time scale, if the star formation process 
took place more efficiently in an enhanced density region, as is the case for n = 2. 

We noticed a difference in the [Fe/H] value at which the abrupt increase of < v° bs t (Z) > 
/cr v (Z) takes place between the Galaxy and M31. By comparing these observed curves with 
the results of simple models, we derived ratios of the collapse time scale to the star formation 
time scale for the Galaxy and M31 of t co i/t s f ~ 3 — 11 and t co i/t s f > 7, respectively, when 
the power index of the Schmidt Law for star formation is n = 1. For n = 2, those ratios are 
tcoi/tsf > 100. This might be a evidence that the n = 2 model is ruled out, however we have 
to construct more realistic chemo-dynamical evolution model for discussing whether it is 
ruled out or not. These differences in the [Fe/H] value do not reflect a significant differences 
in the t co i/t s f. We suggest that the difference of the < v rot > /ovmetallicity relation 
between galaxies should be a probe for studying the differences of formation scenarios of 
those galaxies, although we cannot restrict the formation scenario of a galaxy from it. 

For the GCSs of both our Galaxy and M31, we see that the value of < vfj(Z) > /a v (Z) 
is ~ at —1.5 <[Fe/H]< —0.8, which is not consistent with the models. Since we assumed 
that the protogalactic gas was monolithic in our models, < v° b \{Z) > /a v (Z) may have been 
affected by merger events during the epoch that correspond to that metal range. 

GCs are useful probes for investigating the chemo-dynamic properties of host galaxies. 
New observations of the GCSs of nearby disk galaxies with 8- 10m class telescopes will provide 
fruitful suggestions on the history of galaxy formation. 
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A. Scheme for deriving the rotational velocity < v° b f(Z) > and Velocity 

dispersion a v (Z) 

We describe the scheme we used to derive the average rotational velocity < v° h Q l(Z) > 
and the velocity dispersion a v (Z) from the velocity components (v ^(r, Z), (r,Z), v * x (r, Z))(cf. 
Figure 7) observed line-of-sight velocity v° b s s (r,Z) of GCs. 

The observed line-of-sight velocity vf^ s (r, Z) of GCs is the summation of the projected 
components of the actual 3D velocity components (v^^r, Z), v^(r, Z), v*J-(r, Z))(cf. Figure 
7 for a disk galaxy and Figure 8 for our Galaxy) and as is expressed by 

v? b s s (r,Z) = v t r 1 {r,Z)sm6cosi + vl 1 {r,Z)cos6cosi + v t z 1 {r,Z)smi , (Al) 
for a disk galaxy and 



Vi b o S s(r,Z) = v?(r,Z)Jl - # sin 2 l—^—Wl (r, Z) ^ sin I (r, Z) 



r 2 Vr' 2 + z 2 r Vr' 2 + z 2 z v ' y ^r' 2 + z 2 

(A2) 

for our Galaxy, where the superscript t\ means that the velocity components are those 
observed at present(t = ii)(cf. section 2.1). The average rotational velocity for GCs with 
metallicity Z is defined as 



< v obs 

^ u rot 



(Z)>= [ [ v^ s (r,Z)P(r-9,z,vt\vl\v^)rdrd9dzdv^dv^dv t z \ (A3) 
Jr Jv 

where P(r, 9, z, v I 1 , v 1 ^ , v I 1 ) denotes the probability distribution function of GCs in the phase 
space. We assume, for simplicity, that the probability distribution function P(r,9, z, i** 1 , , v \ x ) 
is spatially uniform and satisfies the symmetry requirements to suppress the net radial or 
vertical motion, i.e., < f* 1 >=< v\ x >= 0. 

Then, the observed rotational velocity is given by 

<v° r b t(Z)>=-<v^(Z)>cost (A4) 
and the velocity dispersion is derived as 

a v (Z) = sj< v? f(Z)>-<v° r b l(Z) > 2 . (A5) 



-14- 



REFERENCES 

Arimoto, N., Yoshii, Y., and Takahara, F., 1992, A&A, 253, 21 
Armandroff, T.E., and Zinn, R., 1988 AJ, 96, 92 
Ashman, K.M., and Zepf, S.E., 1992, ApJ, 384, 50 

Ashman, K.M., and Zepf, S.E., 1998, "Globular Cluster Systems", Cambridge University 
Press 

Barmby, P., Huchra, J. P., Brodie, J. P., Forbes, D.A., Schroder, L.L. and Grillmair, C.J., 
2000, AJ, 119, 727 

Burton, W. B., and Gordon, M. A., 1978, A&A, 63, 7 

Chiba, M., and Beers, T.C., 2000, ApJ, 136, 748 

de Vaucouleurs G., de Vaucouleurs A., Corwin H.G., Buta R.J., Paturel G., Fouque P. 1991, 
" Third Reference Catalogue of Bright Galaxies" , Springer- Verlag, New York (RC3) 

Eggen, O.J., Lynden-Bell, D., and Sandage, A.R. 1962, ApJ, 136, 748 

Fall, S.M., and Rees, M.J., 1985, ApJ, 136, 748 

Frenk, C.S., and White, S.D.M., 1980, MNRAS, 198, 173 

Geogiev, T.B., Bilkina, B.I., and Tikhonov, N.A., 1992, A&AS, 96, 569 

Geogiev, T.B., Tikhonov, N.A., Bilkina, B.I., Getov, R.G. and Nedialkov, P.D., 1993, As- 
trofiz. Issled, 35, 43 

Harris, W.E. 1996, AJ, 112, 1487 (http://physics.mcmaster.ca/Globular.html) 

Hartwick, F.D.A., 1976, ApJ, 209, 418 

Helmi, A., and White, S.D.M., 1999, MNRAS, 307, 495 

Huchra, J. P., Brodie, J.P., Caldwell, N., Christian, C, and Schommer, R. 1996, ApJS, 102, 
29 

Huchra, J.P, Brodie, J.P, and Kent, S.M. 1991, ApJ, 370, 495 
Huchra, J. P., Stauffer, J., and van Speybroeck, L. 1980, ApJ, 259, 57 
Jore, K.P, Broeils, A.H., and Haynes, M.P, 1996, AJ, 112, 438 



- 15 - 

Katz, N., and Gunn, J.E., 1991, ApJ, 377, 365 
Kinman, T.D., 1959, MNRAS, 119, 559 
Kissler-Patig, M., and Gebhardt, K., 1998, AJ, 116, 2237 
Kissler-Patig, M., 2000, astro-ph/0002070 

Maline, G., Hartmann, D.H., and Mathews, G.J., 1991, ApJ, 376, 520 

Mihos,J.G, Walker, I.R., and Hernquist, L., 1995, ApJ, 447, L87 

Minniti, D., 1996, ApJ, 459, 175 

Miyamoto, M., and Nagai, R., 1975, PASJ, 27, 533 

Patsis, PA., and Grosb0l, P., 1996, A&A, 315, 371 

Perelmuter, J., Brodie, J.P, and Huchra, J.P, 1995, AJ, 110, 620 

Perelmuter, J., and Racine, R., 1995, AJ, 109, 1055 

Rees, M.J., and Ostriker, J.P, 1977, MNRAS, 179, 541 

Rubin, V. C, and Ford, W. K., Jr., 1970, ApJ, 159, 379 

Tinsley, B.M., 1980, Fund. Cosm. Phys., 11, 1 

Searle, L., and Zinn, R. 1978, ApJ, 225, 357 

van der Hulst, J.M., 1979, A&A, 75, 97 

Yun, M.S., Ho, P.T.P, and Lo, K.Y., 1993, ApJ, 411, L17 

Zepf, S.E., and Ashman, K.M., 1993, MNRAS, 264, 611 

Zinn, R. 1985, ApJ, 293, 424 

Zinn, R. 1993, in ASP Conf. Ser. 48, The Globular Cluster- Galaxy Connection, ed. J. Brodie 
and G.Smith(San Francisco: ASP), 38 



This preprint was prepared with the AAS IATj^X macros v5.0. 



-16- 



Fig. 1. — Collapse profiles of the scale length for three protogalaxy models. Solid, dashed, 
and dotted lines show the profiles of the exponential (t co ; = lGyr), delayed (t trans = lGyr), 
and early (t trans = O.lGyr) collapse models, respectively. 

Fig. 2. — (a)Evolution of the normalized specific rotation < v° h l(Z) > /o v (Z) as a func- 
tion of the metallicity log(Z/Z Q ) for models with the Schmidt Law Index n = 1 and 
tcoi/tsf = ttrans/tsf = 1- Thin and thick solid lines, which actually overlap, show the expo- 
nential collapse model for t co \ = lGyr and t co \ = lOGyr. Thin and thick dashed lines show 
the delayed collapse model for t tran s = ^Gyr and ttrans = lOGyr. Thin and thick dotted 
lines show the early collapse model for t trans = lGyr and t trans = lOGyr. The line of expo- 
nential collapse with t co i = lGyr overlaps the line of exponential collapse with t co \ = lOGyr. 
(b)Similar diagram, but for t co i/t s f = ttrans/tsf = 10. The lines have the same meaning as 
in Fig. 2(a). 

Fig. 3. — (a)Evolution of the normalized specific rotation < v° h l(Z) > /a v (Z) as a function 
of the increasing metallicity log(Z/Z & ) for models with the Schmidt Law Index n = 2 and 
tcoi/tsf — ttrans/tsf = 1- Other captions are the same as for Fig. 2(a). The line of exponential 
collapse with t co \ = lGyr overlaps with the line of exponential collapse with t co \ = lOGyr. 
(b)Similar diagram, but for t co i/t sf = t tra ns/tsf = 10. 

Fig. 4. — (a)Relation between the observed metallicity index {[Fe/H] = log(Fe/H) — 
log(Fe/ H)q) and the normalized specific rotation(< v° b Q l > /a v ) of the GCS of our Galaxy. 
The observational data are from Harris (1996) (solid circles). The numerical model results 
are as described for Fig. 2(a), but for t co i/t s f = t trans /t s f = 6.0. (b)Similar diagram, but for 
n=2 and t co i/t sf = t trans /t s f = 600. 

Fig. 5. — (a)Relation between the metallicity and normalized specific rotation of the GCS 
of M31. The observational data are from Barmby et al. (2000) (solid circles). The numerical 
model results are as described for Fig. 2(a), but for t co i/t s f = t trans /t s f = 15.0. (b)Similar 
diagram, but for n=2 and t co i/t s f = t tra ns/tsf = 1-5 x 10 4 . 

Fig. 6. — (a)Relation between the metallicity and normalized specific rotation of the GCS 
of M81. The observational data are from Perelmuter, Brodie & Huchra (1995) (solid circles). 
The numerical model results are as described for Fig. 2(a), but for t co \/t s f = tt ra ns/t s f = 15.0. 
(b) Similar diagram, but for n=2 and t co i/t s f = t trans /t s f = 1.5 x 10 4 . 

Fig. 7. — Observed velocity components of GCs of a disk galaxy. The upper figure show the 
disk galaxy seen edge-on. The lower figure show the disk galaxy having inclination < 90°. 

Fig. 8. — Observed velocity components of GCs of our Galaxy. GCs are observed from the 



solar position R = 8kpc of Galactic disk. 
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